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SECTION  1 
INTRODUCTION 


Dynamic  tensile  yielding  and  fracture  have  been  identi¬ 
fied  as  priority  areas  for  materials  research. The  Univer¬ 
sity  of  Dayton  Research  Institute  has  been  active  in  this  field 

(2  3) 

for  several  years.  '  This  interim  report  presents  progress 
in  material  description  and  laboratory  techniques  during  the 
period  of  March  to  September  1981. 

The  University  research  effort  has  focused  on  structural 
metals.  Data  are  reported  herein  on  1020  steel,  RHA  (armor 
plate),  several  alviminum  alloys,  OFHC  copper,  and  nickel  200. 

The  experiments  have  subjected  the  sample  materials  to 
various  forms  of  dynamic  tensile  loading.  Test  geometries 
included  flat  plate  impacts,  expanding  rings,  and  tensile 
Hopkinson  bars.  Some  of  the  flat  plate  experiments  were  modi- 
field  from  those  reported  previously by  the  addition  of  fused 
silica  window  discs.  The  expanding  ring  and  Hopkinson  bar  test 
technique  were  wholly  developed  during  the  present  reporting 
period. 

This  report  presents  tensile  yield  data  for  aluminum, 
mild  steel,  RHA,  and  copper.  Tensile  fracture  data  are  pre¬ 
sented  for  1100  aluminum,  RHA,  and  nickel.  Developments  in 
laboratory  techniques  are  reported  for  every  type  of  test. 


SECTION  2 

SPALL  DAMAGE  IN  1100  ALUMINUM 


2.1  IMPACT  EXPERIMENTS 

Lower  bounds  for  plate  impacts  that  produce  no  damage 
in  1100  aluminum  were  not  found  in  the  previous  effort. 

Two  plate  impact  tests  were  performed  to  obtain  these  limits. 

(3) 

The  plate  impact  technique  used  previously  was 
modified  slightly.  Flyer  plates  had  a  step  on  their  rear 
surfaces,  so  that  each  of  the  two  halves  was  a  different 
thickness.  Nevertheless,  regions  of  one-dimensional  strain 
were  produced  in  the  target  plate,  in  contrast  to  the  wedge 

(4 ) 

flyers  plates  used  by  some  other  investigators. 

A  new  tilt  pin  was  also  designed.  Three  insulated  wires 
were  inserted  through  holes  in  the  target  plates  and  lapped 
flush  with  the  impact  surface.  They  were  positioned  so  as  to 
be  struck  just  inside  the  edge  of  the  flyer  plates.  The  place 
ment  and  closure  of  these  pins  is  considerably  more  precise 
than  the  self-shorting  pins  used  previously.  Tilt  was  less 
than  could  be  resolved  with  the  timing  circuitry:  50  ns.  In 
every  shot,  impact  velocity,  as  measured  between  a  barrel 
shorting  pin  and  target  contact,  was  in  excellent  agreement 
with  streak  camera  records  of  the  projectile  taken  just  prior 
to  impact. 

The  shots  that  were  done,  376  and  383,  employed  d  =  1 
and  2  mm,  and  T  =  4  mm  at  u^  =  211  and  234  m/s,  respectively. 
Little  damage  was  observed.  An  exceptionally  clear  VISAR 
record  was  obtained  from  shot  383,  Figure  1.  Unfortunately, 
the  viewing  area  was  near  the  center,  so  the  release  wave 
was  not  precisely  one-dimensional.  Two  distinct  arrivals  are 
evident.  In  this  shot  the  HEL  is  clearly  defined  at  u  =  20  m/ 
corresponding  to  =  3.6  kbar,  or  Y  =  0.9  kbar. 


Figure  1.  VISAR  Record  From  Shot  383.  Two  quadrature  channels 
Precursor  and  arrival  of  release  wave  are  easily 
distinguished.  Time  marks  are  1  us. 

A  new  technique  was  developed  for  quantifying  the  spall 
damage  in  targets.  A  cross-section  through  the  damaged  region 
was  prepared.  The  distance  was  measured  from  the  center  of 
each  void  to  the  rear  surface  of  the  target  for  a  significant 
fraction  of  the  void-filled  region.  Care  was  taken  to  avoid 
regions  of  the  target  where  the  initial  shock  was  not  planar, 
due  to  the  finite  diameter  of  the  flyer  plate.  The  mean  of 
these  measurements  was  recorded.  Next,  an  incipient  fracture 
line  was  drawn,  through  most  of  the  voids,  representing  our 
judgement  of  the  path  that  the  spall  fracture  would  have 
developed  had  the  fracture  been  complete.  Lastly,  the  percent¬ 
age,  F,  of  the  incipient  fracture  line  actually  occupied  by 
voids  was  measured.  These  data  are  recorded  in  Table  1.  Table 
1  also  includes  data  from  Reference  3.  Several  previous  inves¬ 
tigations  have  defined  spall  separation  as  50  percent  failure 
along  the  spall  plane.  This  effectively  corresponds  to  F  =  50. 


TABLE  1 

SPALL  OBSERVATIONS  IN  1100  ALUMIi<Ux4 


u 

(m/s) 

d 

(mm) 

T 

(mm) 

Shot 

No. 

Description 

x 

s 

(mm) 

B 

211 

1.0 

4.0 

376 

no  damage 

0 

234 

1.0 

4.0 

383 

no  damage 

— 

0 

300 

1.0 

4.0 

272 

connected  voids 

0.83±0.06 

30 

324 

1.0 

4.0 

271 

connected  voids 

NM 

85 

365 

1.0 

4.0 

270 

connected  voids 

0.79±0.05 

70 

472 

1.0 

5.0 

51 

cracks 

0.80±0.11 

75 

590 

1.0 

5.0 

50 

voids 

NM 

NM 

1556 

1.0 

5.0 

52 

voids 

NM 

NM 

211 

2.0 

4.0 

376 

no  damage 

—  — 

0 

234 

2.0 

4.0 

383 

connected  voids 

NM 

'veo 

409 

2.0 

5.0 

47 

voids 

— 

— 

570 

2,0 

5.0 

49 

connected  voids 

1.95±0.09 

30 

894 

2.0 

5.0 

106 

connected  voids 

NM 

NM 

1024 

2.0 

5.0 

48 

connected  voids 

NM 

NM 

There  are  some  anomalies  in  Table  1  in  which  F  is  low  for  high 
velocity  impact;  this  may  be  due  to  void  closure.  In  only  one 
shot,  number  272,  was  the  fracture  less  than  complete.  Many 
authors  have  suggested  the  concept  of  a  critical  void  fracture, 
,  above  which  coalesence  always  occurs,  followed  by  complete 
rupture.  These  results  show  >  30%. 

2 . 2  SWAP  CALCULATIONS  FOR  1100  ALUMINUM 

The  SWAP  code,  as  described  in  Reference  3,  was  modified 

to  include  the  epoxy  material  behind  the  flyer  plate.  The 

effects  of  the  epoxy  were  significant  for  aluminum  flyer 

plates,  although  not  for  others.  The  residual  stress  at  the 

flyer/epoxy  interface  after  release  was  9  kbar.  This  would 

have  been  zero,  had  the  flyer  been  backed  by  air.  A  number 

of  runs  using  various  value  of  Y  (yield  stress)  and  o  (spall 

s 

stress)  were  done  for  shot  272.  Depending  on  the  values  of  Y 
and  Og  input  into  the  code,  spall  took  place  on  either  the 


4 


second  elastic  release  wave  or  the  leading  plastic  release 

wave  reflected  from  the  flyer  rear  surface.  These  cases  are 

illustrated  in  Figures  2  and  3*.  In  Figure  2  spall  occurs 

at  X  =  0.56  nun  and  in  Figure  3  at  x  =1.00  nun.  Actually, 
s  s 

spall  occurs  at  x  =  '\j0.8  mm  (see  Table  1).  Thus,  SWAP 

s 

is  not  capable  of  predicting  the  location  of  the  spall  plane. 

It  seems  unlikely  that  any  reasonable  modification  of  the 
release  wave  profiles  could  affect  this  discrepancy.  The  spall 
signal  must  have  originated  after  void  initiation  because 
it  did  not  originate  at  the  tail  end  of  the  release  fan  from 
the  target  rear  surface.  It  is  along  this  trailing  charac¬ 
teristic  that  any  specified  level  of  tensile  stress  is  first 
achieved.  Thus,  it  appears  that  an  incubation  period  was 

SM0>  r-oirs.PLOts- »L/Ai.o.30on/s  */2i /8i 


Figure  2.  Characteristics  Calculated  by  SWAP  for  Shot  272. 

Y  =  3  kbar,  Og  =  8  kbar.  Impact  is  at  x  =  0, 
t  =  0.  (Flyer  is  backed  by  air.) 


*  In  these  figures,  as  in  other  (x,t)  diagrams  presented  in  this 
report,  impact  occurs  at  this  origin.  The  flyer  plate  moves 
from  left  to  right.  There  is  a  reference  line  near  the  spall 
plane  that  does  not  represent  a  material  interface. 
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Figure  3.  Characteristics  Calculated  by  SWAP  for  Shot  272 
1.0  kbar,  o  =8  kbar. 


necessary  before  the  material  "opened"  on  the  spall  plane, 
generating  the  spall  signal.  Several  different  values  of  Y 
were  tried  in  order  to  match  the  observed  free  surface  velocity 
decay.  Y  =  1  kbar  seemed  to  be  best  which  is  consistent 
with  the  observations  from  shot  383.  A  free  surface  history 
computed  for  shot  272  is  shown  in  Figure  4.  The  free  surface 
history  measured  in  the  VISAR  is  shown  in  Figure  5  (the  time 
scale  has  been  revised  from  Reference  3) .  These  two  velocity 
histories  still  contain  several  descrepancies  that  are  dif¬ 
ficult  to  resolve.  In  the  measurement,  0.25  us  elapsed 
between  the  elastic  precursor  and  shock  arrival.  In  addition, 
the  precursor  was  not  sharply  rising.  The  precursor  is  much 
more  distinct  in  shot  383  (Figure  1) ,  in  which  the  precursor 
arrival  precedes  the  main  shock  by  0.21  us.  If  we  synchronize 
the  VISAR  data  and  SWAP  calculation  by  the  arrival  of  the 
plastic  wave,  then  it  is  clear  that  the  actual  velocity  decay 
is  more  gradual  than  the  predicted  one.  However,  we  are 
not  convinced  that  the  VISAR  data  from  shot  272  is  an 
accurate  portrayal  of  the  free  surface  velocity  variation  in 
an  aluminum  target  undergoing  incipient  spall.  We  intend  to 
perform  several  replications  of  the  shot  before  attempting  to 
reach  final  conclusions.  However,  in  general  analysis  of  the 
discrepancies  between  calculated  and  measured  free  surface 
histories  indicates  that  an  elastic-plastic  model,  as  is  con¬ 
ventionally  used  to  describe  release  waves,  and  as  is  presently 
used  in  SWAP,  is  incapable  of  reproducing  the  observed  velocity 
history  and  spall  plane  location.  The  release  process  is  much 
more  gradual  than  is  predicted  by  the  model,  in  which  two 
rarefaction  shocks  of  strength  Y  propagate  throurh  the  shock 
compressed  material  at  the  speed  of  longitudinal  elastic  waves. 
In  reality,  all  of  the  rarefaction  waves  seem  to  be  dispersive. 

In  the  SWAP  calculations  for  shot  272,  spall  occurred 

at  a  =-12.0  kbar  and  did  not  occur  at  -  -12.5  kbar.  Thus, 
s  s 

the  threshold  spall  stress  for  d  =  1  mm  is  between  -12  and 

(5) 

-12.5  kbar.  Barker  has  successfully  modelled  release  profiles 
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such  as  Figures  2  and  5  using  SWAP.  He  employed  a  yield 
stress  that  decreased  during  the  release  process.  We  have 
not  tried  that  approach  yet. 

In  order  to  correctly  derive  spall  criteria,  the  com¬ 
plete  release  wave  profile  must  be  reasonably  accurately  cal¬ 
culated.  However,  in  experiments  such  as  these,  only  the 
leading  characteristics  of  the  release  fan  can  be  observed, 
because  only  they  arrive  at  the  target  free  surface  undis¬ 
turbed  by  the  spall  process.  Thus,  complete  calibration  of 
a  SWAP  model  is  not  possible  from  these  experiments.  In  order 
to  obtain  better  release  data,  impedance-matching  windows  were 
placed  on  the  target  rear  surface. 

The  revised  target  geometry  is  portrayed  in  Figure  6. 

The  targets  were  backed  with  blocks  of  fused  silica.  Use  of 
such  a  "window"  material  has  been  pioneered  by  researchers  at 


Figure  6.  Target  Design  Using  Fused  Silica  Window.  (Dimensions 
in  mm. ) 
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Sandia  Corporation  (see  Reference  6  for  example) .  Fused 
silica  has  a  shock  impedance  similar  to  aluminxim.  Hence, 
only  a  modest  reflected  tensile  wave  is  produced  when  the 
impact-generated  shock  wave  arrives  at  the  metal/silica  inter¬ 
face.  Spall  does  not  occur  and  the  release  waves  propagating 
from  the  flyer  plate  rear  surface  are  not  significantly  dis¬ 
turbed.  However,  the  VISAR  signal  must  be  corrected  for  the 
effect  of  variations  in  the  index  of  refraction  of  the  silica 
window. 

Reference  6  presents  data  for  the  dependence  of  index  of 
refraction  of  silica  on  stress.  Hence,  in  order  to  reduce  the 
VISAR  data,  a  SWAP  calculation  was  run  to  determine  the  stress 
state  in  the  silica  as  a  function  of  time.  The  compression  in 
the  silica  is  always  elastic,  since  Y  =  98  kbar  in  this  material. 
The  loading  data  of  reference  6  were  used: 

a  =  776e  -  4159e^  +  30340e^  -  69260e^  (1) 

(for  0  in  kbar) .  Equation  (1)  has  the  unusual  property  that 

2  2 

there  is  a  substantial  stress  range  for  which  5  a/de  <  0. 

It  turned  out  to  be  very  difficult  to  run  this  model 
in  SWAP,  because  SWAP  assumes  that  elastic  waves  "shock  up"  and 
propagate  at  Cj^.  The  release  wave  was  also  hard  to  model  be¬ 
cause  9K/9P  <  0  always  led  to  code  failure.  Finally,  it  turned 
out  to  be  practical  to  adjust  the  value  of  the  shear  modulus 
so  that  equation  (1)  was  fulfilled. 

The  results  of  a  sample  calculation  are  shown  in  Figure  7. 
It  is  for  an  impact  velocity  u^  =  418  m/s.  The  peak  stress 
was  33.2  kbar  in  the  aluminum  and  27.5  kbar  in  the  silica.  The 
anomalous  behavior  of  the  sound  speed  in  silica  can  be  seen 
by  examining  the  evolution  of  the  release  waves  emanating  from 
the  rear  surface  of  the  flyer  plate.  They  spread  in  the  alumi¬ 
num,  but  shock  up  in  the  silica. 
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Figure  7.  Characteristics  in  (x,t)  Plane  for  Aluminum  Target 
Backed  by  Si02  Disc. 

Eight  shots  have  been  conducted  with  silica  windows. 
Unfortunately,  none  have  yielded  fully  satisfactory  results. 

The  shots  are  listed  in  Table  2.  No  spall  occurred  in  the 
targets.  In  the  first  four  tests,  the  VISAR  was  run  with 
glass  etalons,  for  which  the  fringe  constant  was  K  =  338  m/s 
fringe. 

The  first  two  tests  were  intended  for  practice  and  used 
nickel  targets.  Polished  silica  windows  were  pressed  directly 
to  the  targets.  The  results  were  encouraging. 

(7) 

Following  practice  at  Sandia  ,  for  subsequent  targets 
the  silica  front  surface  was  depolished  and  aluminized.  It 
was  bonded  to  the  target  with  Hysol  epoxy. 

Of  the  aluminum  targets,  only  shot  385  (Figure  8)  yielded 
good-looking  data.  However,  the  beam  intensity  racord  for  this 
shot  was  lost.  (It  was  off-scale.)  Some  of  the  fluctuations 
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TABLE  2 

SUMMARY  OF  SHOTS  WITH  Si02  WINDOWS 


1 

(m/s) 

d 

(mm) 

'  T 

(mm) 

Material 

Shot 

No. 

Results 

576 

2 

5 

Ni200 

381 

very  noisy 

o 

00 

2 

;  5 

Ni200 

382 

noisy,  but  reducible 

432 

1 

i  4 

llOOAl 

384 

off  scale 

418 

1 

4 

llOOAl 

385 

noisy,  but  reducible 

216 

1 

4 

llOOAl 

483 

HF  noise 

347 

1 

llOOAl 

485 

HF  noise,  BIM  off  scale 

343 

1 

4 

2024A1 

486 

signals  too  small 

of  the  signal  in  Figure  8  may  be  due  to  beam  intensity  varia¬ 
tions,  because  it  seems  impossible  to  account  for  them  with 
a  plausible  interface  velocity  history. 

In  shot  483  the  VISAR  sensitivity  was  increased  to 
K  =  51  m/s  fringe  in  order  to  obtain  more  fringes  and  relieve 
ambiguity  in  the  data  interpretation.  The  data  were  too  noisy 
to  interpret.  It  also  appeared  that  the  data  were  affected 
by  light  emission  from  the  quartz.  In  shot  484  a  gelation 
filter  was  placed  in  the  VISAR.  However,  there  was  still 
excessive  high  amplitude  high  frequency  noise.  In  shot  485 
a  narrow  band  interference  filter  was  used.  The  BIM  still 
went  off  scale.  The  size  of  the  beam  stop  for  the  reflection 
from  the  silica  rear  surface  was  increased  for  shot  486. 
Unfortunately,  the  resultant  signal  was  of  such  low  cimplitude 
that  data  reduction  was  not  possible. 

Experiments  with  silica  windows  will  continue  in  the 
next  reporting  period.  Our  approach  will  be  to  use  higher 
recording  speeds  on  the  oscilloscopes  and  lower  impact 
velocities . 
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Figure  8.  VISAR  Data  from  Silica-backed  Target;  Shot  385.  Time 
marks  are  1  ys.  Two  quadrature  signal  channels  shown. 
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SECTION  3 

SPALL  THRESHOLDS  FOR  RHA 


3 . 1  IMPACT  TESTS 

The  impact  test  data  base  reported  in  Reference  3  was 
expanded  in  order  to  determine  the  dependence  of  spall  stress 
on  tensile  pulse  duration.  The  total  data  set  is  presented 
in  Table  3.  The  stepped  flyer  plates  and  damage  assessment 
techniques  described  in  the  previous  chapter  were  used  for 
many  of  the  new  shots. 

Complete  spall  plate  separation  did  not  occur  in  the 
stepped  flyer  plate  experiments.  This  was  apparently  because 
each  half  of  the  incipient  spall  plate  could  not  separately 
acquire  sufficient  momentum  to  rupture  the  unspalled  material 
around  its  periphery.  However,  in  only  one  case  was  the  damage 
only  partial — indicating  that  coalescence  of  voids  and  cracks 
in  RHA  occurs  very  rapidly  after  their  formation,  and  that 
the  critical  void  fraction  for  rupture  is  relatively  small. 

3.2  SWAP  CALCULATIONS 

The  RHA  impacts  were  modelled  with  the  SWAP  code.  For 
these  calculations,  the  RHA  Hugoniot  from  reference  8  was  used: 

U  =  4.51  +  1.43  u.  (2) 

A  constant  yield  stress  of  4.7  kbar  was  used,  based  on  the 
observed  elastic  precursor  amplitudes. 

For  each  shot,  the  most  tensile  stress,  o  .  ,  was  com- 

min 

puted  at  the  spall  plane  location  (spall  being  suppressed  in 
the  code).  The  results  are  shown  in  Figure  9.  It  can  be 
seen  that  the  data  marginally  support  the  assertion  that  more 
tensile  stress  is  required  for  fracture  of  the  material  when 
the  pulse  length  is  relatively  short.  Spall  may  be  approximately 
characterized  by  a  critical  stress  spall  stress,  o^,  equal 
to  -52  kbar. 
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TABLE  3 

SPALL  THRESHOLD  DATA  FOR  RHA 


(m/s) 

d 

(mm) 

T 

(mm) 

Shot 

No. 

Description 

Xs 

(mm) 

Hi 

391 

0.5 

4.0 

379 

no  damage 

-- 

0 

283 

1.0 

4.0 

371 

no  damage 

— 

0 

319 

1.0 

4.0 

372 

isolated  voids 

1.17+0.13 

1.63+0.13 

55 

355 

1.0 

4.0 

373 

cracks 

1.21+0.01 

95 

391 

1.0 

4.0 

379 

cracks 

1.07+0.03 

90 

235 

2.0 

4.0 

265 

no  damage 

— 

0 

247 

2.0 

4.0 

62 

no  damage 

— 

0 

283 

2.0 

4.0 

371 

no  damage 

0 

319 

2.0 

4.0 

372 

cracks 

2.00+0.14 

1.43+0.15 

90 

355 

2.0 

4.0 

373 

cracks 

2.31+0.03 

95 

,  380 

2.0 

5.0 

60 

spall 

— 

100 

383 

2.0 

4.0 

371 

spall 

100 

390 

2.0 

4.0 

268 

spall 

100 

535 

2.0 

5.0 

25 

spall 

100 

771 

2.0 

4.0 

21 

— 

spall 

mm  mm 

100 

The  Tuler-Butcher  damage  parameter,  ,  was  also  cal¬ 
culated  for  the  RHA  impacts.  The  parameter  is  defined  by 

K,  =  /(o-o  ) ^  dt  (3) 

A  O 

1k^ I  >  Ik^I  for  spall 

The  model  used  to  describe  fracture  in  reference  9  assumed 
that  voids  initiate  at  -25  kbar.  Therefore,  calculations 
were  executed  with  =  -20,  -25,  and  -30  kbar  in  Equation  (3). 
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0.5  I.  1.5  2. 

d(mm) 


Figure  9.  Effect  of  on  d  on  Spall  in  RHA 


0.5  I.  1.5  2. 


d(mm) 

Figure  10.  Effect  of  K2  and  d  on  Spall  in  RHA. 


16 


Both  A  =  1  and  A  =  2  were  tried.  Combinations  were  sought 

that  would  result  in  a  value  of  that  does  not  depend  on 

impact  parameters.  The  best  values  were  Oq  =  -30  kbar  and 

2 

A  =  2,  shown  in  Figure  10.  =  *^120  kbar  -s.  However, 

this  formula  fits  the  data  no  better  than  the  simple  spall 

stress  criterion  with  o  =  -55  kbar.  Therefore,  over  the 

s 

loading  range  examined  here,  there  is  really  no  incentive  to 
use  a  Tuler-Butcher  time  dependent  criterion  for  this 
material . 


1 
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SECTION  4 

SPALL  IN  NICKEL  AND  COPPER 


The  data  base  established  in  Reference  3  for  these 
materials  has  been  expanded.  As  a  consequence,  the  spall 
conditions  are  more  precisely  determined  for  these  materials. 

4.1  SPALL  OBSERVED  IN  NICKEL 

The  material  used  here  was  nickel  200.  This  is  a  com¬ 
mercial  pure  alloy.  The  hardness  of  the  targets  was  in  the 
range  of  67  to  75  on  the  Rockwell  B  scale  (1/lC-inch  diameter 
ball,  100  kg  load) . 

Table  4  shows  the  observations  made  on  impacted  targets. 
The  table  includes  both  new  results  and  those  from  Reference  3. 
Figure  11  is  a  plot  of  maximum  compressive  stress  (very  nearly 
equal  to  a  .  )  versus  flyer  plate  thickness  for  Ni-200.  SWAP 
calculations  for  this  material  have  not  been  performed  yet. 

The  data,  as  plotted  in  Figure  11,  indicate  a  strong 
dependence  of  spall  stress  on  load  duration.  However,  these 
results  must  be  interpreted,  with  some  caution  because  two 
different  batches  of  nickel  were  used.  All  the  specimens 
except  shots  374  and  375  were  from  one  batch,  for  which  Figure 
12  shows  sample  grain  structure.  Data  from  these  two  shots 
are  flagged  in  Figure  11.  Shots  274  and  275  used  material 
(also  supplied  and  certified  as  Ni-200)  that  had  a  slightly 
smaller  grain  structure  (see  Figure  13) ,  although  the  same 
hardness,  as  the  other  material.  However,  the  274  and  275 
material  also  contained  inclusions  (also  shown  in  Figure  13) . 
Failure  seemed  to  initiate  preferentially  at  the  inclusions. 

In  shot  274,  although  there  was  no  preponderance  of  damage 
on  the  spall  plane,  there  were  several  voids  that  had  clearly 
initiated  at  inclusion  sites.  Consequently,  the  included 
material  may  be  supposed  to  spall  more  easily.  This  explains 
why  there  is  a  more  damage  in  shot  375  than  shot  119. 
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Consistent  with  this  interpretation,  we  can  confidently  retain 
shot  274  as  indicating  impact  conditions  that  do  not  result 
in  spall  in  nickel  200. 

Spall  from  thick  flyer  plates  (d  _  '2  mm)  can  be  pre¬ 

dicted  by  setting  =  -30  kbar.  However,  spall  induced  by 
thinner  flyer  plates  will  require  a  history  dependent  fracture 
model . 

In  shot  382  free  surface  release  was  measured  for  a  nickel 
target  backed  by  a  fused  quartz  v;indow  (see  Section  2)  .  The 
VISAR  trace  was  noisy,  but  an  interpretation  seemed  feasible. 

An  HEL  of  5.6  kbar  is  observed,  corresponding  to  Y  =  6.6  kbar. 
This  is  consistent  with  Hopkinson  bar  measurements  of  tensile 
strength.  On  the  one  specimen  tested,  Y  increased  from  an 
initial  value  of  2.5  kbar  to  8.5  kbar  at  20  percent  strain. 


Figure  13.  Microstructure  of  Ni-200  Used  in  Shots  274  and  275, 
Showing  Inclusion. 


4.2  RESULTS  FOR  COPPER 

The  data  reported  in  reference  3  showed  the  impact  con¬ 
ditions  resulting  in  incipient  spall  fracture  in  OFHC  copper 
(approximately  half-hard).  However,  there  was  an  anomaly 
in  the  data,  in  that  one  shot  (number  109)  utilized  an 
aluminum  flyer  plate.  The  results  of  that  shot  were  not  con¬ 
sistent  with  the  other  data  obtained  with  copper  flyer  plates. 

In  order  to  resolve  this  conflict,  shot  377  was  conducted, 
utilizing  an  1100  aluminum  flyer  plate  with  1.5  mm  and  2.0  mm 
steps.  The  results  are  shown  in  Table  5,  together  with  the 
data  of  reference  3.  The  new  results  are  much  more  consistent 
than  shot  109.  We  conclude  that  shot  109  is  fallacious;  per¬ 
haps  the  velocity  was  incorrectly  measured. 

The  most  negative  stress  in  shot  377  was  -41  )<bar, 

which  is  less  than  o  ,  according  to  reference  3  (Figure  20). 

s 

The  damage  consisted  of  many  closely  spaced  voids  that  would 

presumably  would  have  formed  a  crac)c  had  a  constant  thic)cness 

flyer  plate  been  used.  Figure  14  shows  K2  for  the  data  in 

Table  5  (without  shot  109)  .  It  is  evident  that  these  results 

2 

are  consistent  with  50  <  <  110  )cbar  -  ijs,  as  concluded 

in  reference  3.  (Fee  reference  3  for  a  discussion  of  the 
datum  at  d  =  0.4  mm.) 
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TABLE  5 

SPALL  IN  COPPER 


1 

Velocity 

(m/s)  1 

1 

Max  i 
Stress  ! 
(kbar)  j 

Flyer 

Thickness 

(mm) 

Target 

Thickness 

(mm) 

Shot 

No. 

Remarks 

i 

80  I 

14 

1.50 

4.00 

121 

IJo  voids 

103  ! 

i 

19 

1.45 

3.90 

111 

Void  layer 

244 

45 

1.45 

4.00 

110 

Spall 

separation 

366 

41 

1.50  A1 

4.00 

377 

Void  layer 

168 

30 

2.0 

4.00 

221 

Spall 

separation 

189 

35 

2.0 

4.00 

220 

Spall 

separation 

322 

35 

3.0  (Al) 

3.96 

109 

Void  layer 

366 

41 

2.0  (Al) 

4.00 

377 

Void  layer 

286 

53 

1 

2.5 

3.94 

108 

Spall 

separation 

Spall 

separation 


334 


62 


2.5 


4.0 
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SECTION  5 

EXPANDING  RING  TESTS 


5 . 1  BACKGROUND 

Only  a  few  test  procedures  are  available  for  determination 

of  constitutive  relations  at  negative  mean  stress  and  high 

4  -1 

strain  rates.  The  range  around  10  s  is  particularly  dif¬ 
ficult  to  achieve.  Rapid  mechanical  machines  are  unable  to 

2  -1 

achieve  tensile  stresses  above  10  s  .  The  split  Hopkinson 
bar  tensile  test  (see  next  section)  can  extend  this  range 
another  decade.  Resolution  of  rapidly-changing  stresses  is 
difficult  with  a  Hopkinson  bar  because  of  the  effects  of  the 
sample  grip  on  the  transmitted  waveforms. 

Most  studies  of  tensile  failure  at  high  strain  rates 
utilize  i.vceraction  of  tensile  waves  produced  from  impact¬ 
generated  shock  waves,  as  in  the  spall  work  reported  in  pre¬ 
ceding  sections.  These  techniques  are  limited  in  that  the 
subject  material  experiences  intense  compressive  loading  prior 
to  tensile  loading,  and  only  states  of  one-dimensional  strain 
are  achieved. 

The  expanding  ring  technique  presents  an  attractive 
supplement  to  conventional  tensile  test  techniques.  Relatively 

4  _  1 

high  rates  can  be  achieved,  >10  s  .  Data  analysis  is  free 
from  complications  involving  wave  propagation.  Eventually, 
triaxial  data  may  be  obtained  by  testing  rings  with  pre¬ 
drilled  holes. 

In  the  expanding  ring  test,  an  axially-symmetric  thin 
short  ring  of  material  is  driven  radially  outward.  The  ring 
may  be  launched  explosively ,  with  a  magnetic  field^^^^, 
with  an  exploding  wire^^^^,  or  via  a  driver  cylinder . 

Most  past  work  with  expanding  rings  has  been  designed  to 
measure  ultimate  strain  at  fracture.  The  ultimate  strain  is 
given  by 
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Where  the  average  velocity,  t,  is  monitored,  strain  rate  may 
be  calculated  from 


t  =  f/r  (5) 

Here  r^  and  r^  are,  respectively,  the  initial  radius  of  the 
cylinder  and  the  radius  at  failure. 

In  the  experiments  of  reference  15  as  well  as  in  those 
developed  by  the  University,  ring  acceleration  can  be  precisely 
measured.  Before  considering  the  expermental  technique  for 
generating  the  expansion,  it  is  useful  to  consider  the  relation 
ship  between  stress,  strain,  and  ring  motion.  Analytical 
equations  are  evaluated  with  sample  numerical  values  correspond 
ing  to  a  particular  system  of  interest  -  a  20  mm  diameter,  1  mm 
thick  copper  ring  (density,  p^,  of  8.82  g/cm^) . 

The  equations  of  equilibrium*  for  the  expanding  ring  are: 


(6) 


(7) 


1 

r 


0 


and  the  boundary  conditions  are; 


o  =  0  at  outside  surface 
r 

=  P  at  inside  surface 

o,  0 
z 


(8) 


and  are  respectively,  stress  in  the  radial,  azi¬ 
muthal  and  axial  directions.  P  is  the  internal  driving 
pressure . 


In  order  to  calculate  from  f,  it  is  necessary  for  P 
to  be  negligible.  In  that  case,  =  0,  and  a  condition  of 
uniaxial  stress  exists.  The  term  9a^/9r  is  of  order  P/h, 
where  h  is  the  instantaneous  wall  thickness.  Since  h  <<  r, 
consideration  of  the  first  equilibrium  equation  leads  to  the 
requirement 


The  above  inequality  imposes  severe  limitations  on  P. 

We  expect  2  <  Og  <  4  kbar.  Take  3  kbar  for  illustrative 
purposes.  The  ratio  h/r  will  be  typically  between  0.1  and 
0.01.  Lets  take  0.05  as  typical.  Then  we  require  that  P  <<  150 
bars.  This  condition  is  very  difficult  to  achieve  with  explo¬ 
sives;  the  pressure  of  an  expanding  explosive  may  be  estimated 
from; 


P  (in  bars)  =  (^('' 

o 

where  y  is  the  ratio  of  specific  heats  and  is  a  property 

of  the  explosive  and  is  equal  to  the  volume  which  the  detonation 
reaches  in  an  adiabatic  expansion  to  1  bar.  is  the  STP 

specific  volume  of  the  detonation  products.  Typical  values  of 

3 

are  'V'600  cm  /g,  and  a  typical  value  for  y  at  low  temperatures 
is  1.4.  If  the  strain  is  100  percent,  and  blowby  is  ignored, 
it  turns  out  that  P  =  375  bars.  This  is  several  times  higher 
than  can  be  tolerated  by  the  assumption  is  negligible. 

Thus,  direct  explosive  driving  of  rings  and  cylinders  is  not 
an  adequate  technique  for  present  purposes . 

Most  of  the  energy  in  the  cylinder  is  kinetic.  If  a 
strain  rate  is  to  be  obtained  at  a  particular  radius  r^, 
then  the  kinetic  energy  per  unit  length,  T' ,  is: 


T' 


(11) 
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For  the  ring  in  question  at  f  =  5  x  10  s  and  e .  =  In  r./r  = 

10  J  J  ° 

1,  this  equation  gives  1.39  x  10  ergs/cm.  The  strain  energy 

in  the  cylinder  at  failure  is  approximately  Ye^,  where  Y  =  the 

9  ^3 

yield  stress.  This  is  equal  to  3.9  x  10  ergs/cm  for  = 

9  ^ 

1.3  and  Y  =  0.3  kbar,  or  1.23  x  10  ergs/cm  of  length.  This 

energy  is  less  than  10  percent  of  the  kinetic  energy.  The 

important  conclusions  here  are:  (1)  in  order  to  achieve  a  given 

strain  rate,  the  cylinder  must  be  given  a  critical  energy  in 

4  -1 

an  impulsive  way;  and  (2)  for  strain  rates  of  >10  s  ,  failure 
will  occur  with  a  minimal  decrease  in  kinetic  energy. 

The  strain  in  the  ring  can  be  determined  from  the  instan¬ 
taneous  diameter.  The  condition  P  %  0  permits  the  stress  in 
the  ring,  Og,  to  be  determined  by  measuring  the  radial 
deceleration. 


Og  = -prf  (12) 

If  r  is  determined  over  intervals  as  Ar/At,  then  the 
uncertainty  in  r,  5r,  may  be  expressed  as 

6r  =  {(5Af-/At)^  +  (Af6At/At^)  (13) 

As  a  numerical  example,  consider  o„  =  3  kbar  and  r  =  1.25  x 

8  2  ^ 

10  cm/s  .  If  £■  is  equal  to  one-half  a  fringe  at  maximum 

VISAR  sensitivity,  then  Af  =  0.01  ram/us.  The  resolution  should 
be  1/10  of  this,  so  6Af/Af  =  0.1.  This  occurs  over  At  =  8  ys. 
The  uncertainty  in  the  time  interval  will  be  negligible.  Thus, 
6r/r  10  percent.  Since  there  are  two  VISAR  signals  in  quadra¬ 
ture,  the  net  uncertainty  will  be  close  to  7  percent.  The  error 

in  r  is  less  than  '^-0.02  mm  out  of  '^2.7  mm.  Thus,  the  hoop 

stress,  Og,  will  be  determined  to  better  than  10  percent  over 
'V'lO  us  intervals.  The  uncertainty  will  be  less  for  longer 
times  due  to  signal  averaging. 
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5.2  EXPERIMENTAL  TECHNIQUE 

The  technique  used  by  the  University  for  expanding  rings 
uses  an  explosive  charge  inside  of  a  cylindrical  driver  cylinder. 
After  researching  the  availability  of  small  explosive  charges, 
it  was  found  that  significant  cost  and  time  savings  could  be 
achieved  by  using  small  charges  already  on  hand.  These  were 
lead  azide  spheres,  10.8  mm  diameter,  of  1.7  g  mass.  They 
are  centrally  detonated  by  an  embedded  bridgewire. 

An  analysis  was  conducted  to  support  the  design  of  the 
expanding  ring  test.  Equations  were  derived  for  the  expansion 
of  a  cylinder  driven  by  isentropic  expansion  of  an  explosively- 
generated  gas.  If  the  cylinder  inner  radius  is  r^,  and  the 
charge  is  confined,  the  gas  pressure  is: 

P  = 

■^o 

where  P^  is  the  initial  Chapman-Jouget  pressure.  Differential 
equations  for  cylinder  motion  were  derived  and  numerically 
solved.  A  sample  solution  is  shown  in  Figure  15,  in  which 
calculations  are  compared  with  data  from  reference  16.  At 
least  at  early  times,  the  agreement  is  satisfactory. 

Seventeen  explosive  tests  were  conducted  in  order  to 
develop  the  expanding  ring  technique.  The  first  seven  tests 
were  dedicated  to  developing  a  driver  cylinder  that  would 
undergo  large  radial  expansion  without  fragmentation  or  gas 
leakage.  The  above  analysis  guided  the  design.  Three 
materials  were  used;  low  density  polyethylene,  mild  steel, 
and  6061-T6  aluminum. 

The  first  two  tests  were  conducted  with  polyethylene  to 
determine  the  effect  of  cylinder  length.  The  cylinders  were 
64  mm  OD  and  11.9  mm  ID.  50  mm  and  25  mm  length  tubes  resulted 
in  identical  expansion,  viz.  0.36  mm.  All  subsequent  cylinders 
had  a  length  of  50  mm.  Test  6  utilized  a  smaller  diameter 
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RADIAL  EXPANSION  R-R^  (mm) 


12 


Figure  15.  Expansion  of  Explosive  Filled  25  nun  Diameter 
Copper  Cylinder. 
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polyethylene  cylinder  in  order  to  increase  the  expansion.  The 
OD  was  50  nun.  The  radial  expansion  more  than  doubled,  but 
was  still  only  0.86  mm. 

Tests  3,  4,  and  5  were  with  steel  cylinders.  The  ID's 
were  11.9  mm.  The  variation  of  expansion  with  OD  can  be  seen 
in  Figure  16.  The  rapid  increase  in  Ar/r  with  decrease  r  led 
us  to  believe  that  cylinder  significantly  smaller  than  18  mm 
diameter  would  rupture. 

An  aluminum  cylinder  was  used  in  Test  7.  It  had  an  OD  of 
25.5  mm.  Although  the  cylinder  remained  intact,  radial  cracks 
were  produced,  and  it  was  clear  that  even  a  very  small  reduction 
in  diameter  would  result  in  shatter. 

The  next  stage  in  development  of  a  ring  test  was  to 
demonstrate  that  rings  could  be  launched  in  free  expansion  at 
velocities  of  interest.  For  the  developmental  shots,  it  was 
decided  that  a  streak  camera  would  be  adequate  for  velocity 
measurements.  The  camera  employed  was  a  Gordin  Model  231, 
which  uses  70  mm  wide  film  and  is  continuous  access. 

o  STEEL,  NO  RtNG 


Figure  16.  Relative  Growth  in  Driving  Ring  Diameters.  Cross 
indicates  partial  fracture. 
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Four  shots  were  done  to  develop  experimental  techniques. 

The  most  difficult  aspect  was  to  assure  prompt  charge  initiation, 
consistent  with  the  limited  recording  time  of  the  optical  data 
train . 

The  system  eventually  developed  is  shown  in  Figure  17.  A 
Vivitar  flash  unit  was  used  for  illumination.  The  flash  dura¬ 
tion  was  adjusted  to  be  600  ;is.  The  flash  had  an  unavoidable 
200  ys  pause  in  its  trigger  circuit.  The  camera  was  adjusted 
to  run  at  between  0.5  and  1  mm/ys,  so  that  the  entire  flash 
duration  could  be  recorded  without  rewrite.  The  explosive  had 
an  unpredictable  jitcer  of  several  hundred  microseconds,  even 
with  steeply-rising  high  voltage  initiation  pulses.  This 
jitter  resulted  in  occasional  loss  of  data.  Figure  18  shows 
a  sample  streak  record  of  a  detonating  charge.  The  initial 
expansion  velocity  is  almost  7  mm/ys. 

In  shots  12  and  14  steel  cylinders  were  employed  to 
launch  copper  rings.  The  rings  were  cut  from  as-received 
OFHC  copper  tubing.  In  Shot  12,  the  driver  cylinder  OD  was 
18,0  mm  and  ID  wa.s  11.7  mm.  A  1 . 9  mm  thick  ring  was  placed 
on  the  cylinder.  The  ring  expanded  in  radius  0.62  mm,  and 
the  average  logarithmic  strain,  from 

t  =  In  JT^/r^  (15) 

was  0.055  +  0.005.  For  Shot  14  the  ring  thickness  was  reduced 
to  1.0  mm.  However,  the  steel  cylinder  was  a  little  larger, 

OD  =  20.0  mm.  The  ring  only  expanded  0.45  mm,  and  the 
logarithmic  strain  was  0.040  +  0.005. 

In  Shots  15  and  17  aluminum  driver  cylinders  were 
evaluated.  The  shots  were  virtually  identical;  in  Shot  15  the 
OD  was  29.4  mm,  and  in  Shot  17  it  was  30.5  mm.  Copper  rings 
1.96  mm  thick  were  used.  The  average  strains  in  these  tests 
were  0.125  and  0.101,  more  than  twice  the  best  value  obtained 
with  steel  cylinders.  The  aluminum  was  indented  about  0.5  mm 
beneath  the  rings,  which  suggest  that  the  acceleration  was 
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Figure  17.  Layout  for  Streak  Camera  Photographs  of  Expanding 
Rings . 


Figure  18.  Streak  Record  of  Detonating  Charge. 
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gradual  (since  there  was  time  for  plastic  flow  around  the  ring) . 
Unfortunately,  due  to  jitter  in  explosive  detonation  time,  ring 
velocity  data  were  not  obtained  in  either  test.  A  curious 
aspect  of  the  tests  was  that  the  expansion  of  the  aluminum 
was  greater  than  in  Shot  7,  even  though  there  was  no  apparent 
radial  fracture.  This  is  shown  in  Figure  16. 

The  final  test  was  with  a  polyethylene  driver  cylinder. 

The  polyethylene  OD  was  50.5  mm.  An  alviminum  ring  was  used  of 
unknown  composition;  its  thickness  was  3.0  mm.  The  average 
strain  was  only  0.055  +  0.003.  Figure  19  shows  the  cylinder 
after  the  test.  It  is  apparent  that  the  ring  would  have  been 
much  more  effectively  launched  were  it  thinner.  The  plastic 
ruptured  in  one  location  behind  the  ring,  releasing  explosive 
gases  that  would  have  invalidated  a  measurement  of  deceleration. 

5.3  FURTHER  DISCUSSION  OF  ALUMINUM  DRIVER  CYLINDERS 

The  aluminum  driver  cylinders  used  in  Tests  15  and  16 

appear  to  be  strained  to  the  fracture  limit,  based  on  the 

following  analysis.  If  =  0  and  there  is  no  net  volume  change, 

-0  =  Near  the  surface,  and  probably  through  the  cylinder, 

Eg  =  ^/r,  where  (5  is  the  displacement.  The  plastic  strain  is 

e  =  ^0  =  6/r.  For  Test  15,  this  gives  e  =0.14. 

The  handbook  value  for  maximum  elongation  for  6061-T6  aluminum 

is  0.12.  In  a  tensile  test,  e  =  e  • 

z  p 

The  maximum  allowable  e  in  the  cylinders  will  be  less 

p 

than  in  a  tensile  test  because  P/Y  is  greater  in  the  cylinder 

test.  To  see  this,  consider  that  in  the  cylinder  a  =  'vO  and 

o  =  ao,,  where  a  <  1.  Hence,  for  the  cylinder, 
z  0 

P/Y  =  -t=====  (16) 

/2[l+2(l+a) 

For  a  <  1  this  is  greater  than  1/3,  which  is  the  value  of  P/Y 
in  a  tensile  test.  Strain  rate  effects  in  this  material  are 
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Figure  19.  Left:  Top  View  of  Polyethylene  Driven  Cylinder 
and  Expanded  Aluminum  Ring.  Right:  Side  View. 

not  significant  below  about  2,000  s  and  it  is  unlikely  that 
strain  rates  are  higher  than  that  while  plastic  flow  is  taking 
place.  Hence,  it  follows  that  the  aluminum  cylinders  in  Tests 
15  and  17  could  not  have  undergone  significantly  more  expansion 
without  rupture. 

Unfortunately,  the  launch  velocity  for  the  copper  rings 
driven  by  aluminum  cylinders  was  not  measured  in  these  tests. 
However,  it  may  be  estimated  from  the  final  diameter  of  the 
rings  and  the  assumption  that  the  initial  kinetic  energy  of  the 
ring  was  transformed  into  plastic  work. 

The  plastic  work  per  unit  volume  is  approximately  equal 

2 

to  1/2  f'pY*  The  kinetic  energy  per  unit  volume  is  1/2  pu^  , 
where  u^  is  the  launch  velocity.  Thus, 

u  / c  Y / p  (17) 

o  -  p 
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Substitution  with  Y  =  2  kbar  gives  u^  100  ra/s.  The  uncertainty 
in  Y  and  the  plastic  work  probably  yields  a  50  percent  uncer¬ 
tainty  in  u^.  The  strain  rate  is  given  by 

e  =  u  /r  (18) 

o 

4  -1 

Thus,  t  is  in  the  range  0.4  to  1.2  x  10  s  .  The  upper  end 
of  this  range  is  fully  adequate  for  determination  of  useful 
yield  strength  data.  However,  larger  values  of  £  are  required 
if  ultimate  strain  is  to  be  determined. 

Thinner  copper  rings  may  be  launched  at  slightly  higher 
velocities.  If  the  time  of  arrest  of  the  driver  expansion  was 
the  same  under  the  ring  and  away  from  the  ring,  then  the  poten¬ 
tial  gain  in  velocity  by  using  thinner  rings  is  simply  the 
ratio  of  Z\r  alongside  the  ring  and  Ar  directly  under  the  ring. 

This  ratio  turns  out  to  be  1.05.  Thus,  a  very  thin  ring  would 
be  launched  at  about  5  percent  higher  velocity  than  the  ring 
actually  employed  in  Tests  15  and  17. 

5 . 4  PLANS  FOR  FUTURE  TESTS 

Experiments  were  halted  after  Test  17  so  that  the  results 
could  be  assessed  and  new  cylinders  fabricated.  The  following 
conclusions  have  been  inferred  from  the  results. 

•  Lower  impedance  drivers  like  aluminum  and  polyethylene 
are  better  than  steel.  They  provide  a  more  gentle 
launch  of  the  ring  because  of  plastic  flow,  and  they 
are  capable  of  imparting  larger  final  strains  to  the 
ring.  The  preshocking  of  the  ring  noted  by  some  other 
investigators  may  be  avoided. 

•  The  ring  and  driver  configurations  used  in  Test  15 
and  17  are  probably  near  optimum.  The  launch  velocity 
should  be  determined  to  verify  that  the  strain  rate  is 
at  least  5,000  s”^. 
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•  More  experiments  should  be  done  with  polyethylene 

drivers  using  thinner  rings.  The  objective  should  be 
to  launch  the  rings  at  higher  velocities. 
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SECTION  6 


HOPKINSON  BAR  APPARATUS  AND 
RESULTS  OF  MEASUREMENTS 

6.1  HOPKINSON  BAR  APPARATUS 

The  Hopkinson  Bar  apparatus  used  in  this  program  is 
essentially  identical  to  that  described  in  reference  17.  A 
number  of  minor  mechanical  modifications  were  carried  out.  The 
bars  were  realigned,  which  necessitated  reworking  the  teflon 
bushing  mounts,  an  oven  was  designed  and  fabricated,  and  new 
strain  gauges  were  mounted.  A  high  speed  camera  set-up  was 
designed  and  deployed  to  photograph  notched-tensile  specimens 
under  test. 

The  most  substantial  modifications  of  the  bar  from  that 
described  in  reference  17  are  in  the  way  that  data  are  recorded. 
Strain  data  from  each  strain  gauge  are  now  recorded  directly  on 
a  digital  oscilloscope.  The  sampling  rate  is  usually  1  MHz. 

The  gain  is  adjusted  so  that  the  signals  of  interest  are  cap¬ 
tured  with  reasonable  precision,  even  when  this  means  that 
preceding  wave  arrivals  are  off  scale.  A  computer  program 
was  written  to  analyze  the  sti;iin  gauge  records.  The  program 
calculates  strain,  stress,  strain  rate,  and  stress  versus 
strain. 

The  essential  equations  used  to  calculate  specimen  stress 
and  strain  have  been  described  by  many  authors  (for  example, 
reference  17) .  A  brief  account  is  given  here. 

'  Strain  gauges  are  placed  equi-distant  from  the  sample. 

The  striker  bar  is  launched  against  the  pressure  bar,  producing 
a  travelling  compressive  wave.  In  the  present  system,  the 
duration  is  300  ys.  When  the  bar  apparatus  is  used  as  a 
compression  test,  the  initial  compressive  pulse  is  transmitted 
to  the  sample.  The  incoming  pulse  is  partially  transmitted 
through  the  sample,  according  to  the  relationship; 
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(19) 


sp 


u  A,  Z,  Z 
o  b  b  s 

A^Z,  +A  Z 
b  b  s  s 


where  is  the  striker  impact  velocity,  Aj^  and  A^  are  the  bar 
and  sample  cross  section  areas,  and  and  are  the  bar  and 
striker  acoustic  impedance  (e.g.,  Oj^  =  1/2  u^Zj^)  .  The  energy 
that  is  not  transmitted  to  the  specimen  is  reflected.  Thus, 
very  compliant  or  small  samples  (A  Z  small)  transmit  little 
stress,  and  vice  versa.  The  pressure  bar/specimen  interface 
acquires  a  velocity  u  (u  =  o  /Z  )  which  essentially  imposes 

S  3  S  p  S 

a  boundary  condition  of  approximately  constant  strain  rate  on 
the  specimen.  The  stress  in  the  specimen,  is  transmitted 

to  the  second,  or  transmission  bar.  The  particle  velocity 
imparted  to  the  transmission  bar  is: 


2  u  A  Z 

u  =  _ S-S 

^t  A_Z_+A^Z 


S  S 


t  t 


(20) 


Since  A  >  A  ,  and  usually  Z  >  Z  ,  this  transmission  is  often 
quite  weak.  The  reflected  wave  bounces  back  and  forth  in  the 
sample,  diminishing  at  each  reflection  from  the  transmission 
or  pressure  bar.  The  period  is  typically  about  8  us,  so  about 
16  tis  is  required  for  the  specimen  to  come  to  stress  equili- 
brixom  with  the  Hopkinson  bars. 

When  the  bars  are  used  for  tensile  tests,  a  collar  is 
placed  around  the  specimen.  The  collar  supports  the  initial 
compressive  stress.  The  transmitted  wave  reflects  as  a  tensile 
wave  at  the  end  of  the  second  bar.  When  the  tensile  wave 
returns  to  the  specimen,  the  collar  falls  away,  and  the  stress 
is  sustained  by  the  specimen.  Figure  20  shows  an  (x,t)  wave 
diagram  for  a  tensile  test.  It  is  essential  that  the  strain 
gauges  be  positioned  so  that  the  stress  waves  reflected  and 
transmitted  from  the  sample  are  recorded  without  interference 
from  other  waves. 
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Figure  20.  (x,t)  Diagram  for  Tensile  Hopkinson  Bar  Test  (from 

Reference  17) . 

The  assumption  that  stress  in  the  sample  is  uniform 
leads  to  the  equation: 

£g  =  ^  /  tr  dt  (21 

where  c  is  the  wave  speed  in  the  bar,  I  is  sample  length,  and 

e  and  e  are  reflected  strain  and  specimen  strain.  (Here, 
r  s 

as  elsewhere  in  this  report,  stress  and  strain  are  positive 
in  compression.)  The  specimen  stress  is  given  by; 


0  .  E  c 

sp  t 


when  E  and  are  the  Young's  modulus  of  the  bar  and  the  trans 
mitted  strain.  The  specimen  strain  rate  is  given  by; 


t 


2c 

i. 


(23) 


r 


These  equations  show  that  the  specimen  strain  is  obtained  by 
integrating  the  strain  pulse  reflected  from  the  specimen,  the 
specimen  stress  is  proportional  to  the  transmitted  strain  pulse, 
and  the  strain  rate  is  proportional  to  the  reflected  strain 
pulse.  The  equations  show  that  high  strain  rates  are  obtained 
by  high  striker  velocities  or  short  sample  lengths.  The  equa¬ 
tions  are  not  valid  for  elastic  response  because  equilibrium  is 
not  established  during  the  time  needed  for  the  sample  to  "ring 
up"  to  the  bar  boundary  conditions,  as  discussed  above. 

Equations  (21)  and  (22)  give  average,  or  engineering, 
stress  and  strain.  True  stress  and  strain  can  be  computed  from: 


=  o  (1  - 
sp 


(24) 


and 

=  «-n(l  -  e^)  (25) 

The  bar  strain  gauges  are  easily  calibrated  by  performing 
a  compressive  test  with  no  specimen.  Under  those  conditions, 
the  strain  in  the  bars  is  just 

=  uy2c  (26) 

Since  the  initial  loading  pulse  is  of  duration  x  (deter¬ 
mined  by  the  length  and  wave  speed  of  the  strike  bar) ,  the 
maximum  strain  that  can  be  produced  is 


e 


s  ,max 


tx 


(27) 


For  the  bar  used  here,  x  =  300  ys. 

The  maximiun  allowable  striker  bar  velocity  is  the 
velocity  which  causes  the  bar  to  deform  plastically.  For  most 
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tests  the  bars  were  1.27  mm  diameter  4340  steel  at  RC  47.  They 

have  a  yield  strength  of  Y  =  14  kbar,  which  corresponds  to  a 

maximum  strain  of  ■  ,  =  0.007.  It  follows  that  u 

b,max  o,max 

71  m/s.  In  practice,  the  maximum  velocity  that  can  be  achioveci 
is  limited  by  the  dimensions  of  the  striker  bar  and  striker 
bar  guide  system  to  about  25  m/s.  This  means  that  for  a  stiff 
sample,  say  of  4340  steel,  the  maximum  strain  rate  is  about 
600  s  In  a  very  compliant  sample,  the  maximum  strain  rate 

is  2800  s 

Figure  21  shows  a  photograph  of  the  Hopkinson  bar  as  set 
up  in  the  laboratory.  Data  that  have  been  obtained  with  the  bar 
are  presented  and  discussed  in  the  following  sections.  Results 
have  been  very  satisfactory  with  one  exception.  At  higher 


Figure  21.  Photograph  of  Hopkinson  Bar  Apparatus. 
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launch  velocities,  the  strain  pulses  are  noisy.  This  becomes 
a  problem  for  high  strain  rate  studies  (>1000  s  of  low 
impedance  specimens  (such  as  aluminum) ,  for  which  the  trans¬ 
mitted  pulse  is  quite  feeble  and  degraded  by  the  bar  noise. 

The  noise  is  believed  to  be  caused  by  flexural  vibrations  in 
the  striker  bar  and  ringing  in  the  strain  gauge  amplifiers. 

6 . 2  RESULTS  OF  HOPKINSON  BAR  TESTS 

A  summary  of  the  tensile  and  compressive  tests  conducted 
are  presented  in  Table  6.  Seven  different  materials  were 
tested  in  tension  and  one  in  compression. 

6.2.1  Mild  Steel 

The  1020  steel  samples  were  machined  at  the 
University  and  were  tested  in  tension  at  strain  rates  ranging 
from  350  to  1,000  s  A  typical  stress-strain  plot  from  the 

digital  data  is  presented  in  Figure  22.  The  stress-strain 
plots,  for  the  various  strain  rates  tested,  were  smoothed  and 
are  presented  in  Figure  23.  The  dotted  lino  in  Figure  22 
represents  a  smoothed  curve  drawn  to  fit  the  data.  Figure  23 
clearly  indicates  a  strain  rate  dependency  for  the  1020  steel. 
The  yield  stress  increases  with  strain  rate.  Similar  obser¬ 
vations  were  reported  in  reference  18  for  a  mild  steel  tested 
in  tension  in  a  flywheel  impact  machine  at  70°F.  The  "over¬ 
shoot"  often  observed  in  mild  steel  at  the  beginning  of  the 
stress-strain  curve  (reference  17)  could  not  be  distinguished 
from  the  digital  stress-strain  plots,  possibly  because  of  the 
superimposed  oscillations  at  the  beginning  of  the  signal  attri¬ 
buted  to  the  wave  propagation  mechanism  in  a  long  bar.  For  the 
same  reason,  an  accurate  yield  stress  can  not  be  measured;  an 
average  yield  stress  was  estimated  and  reported  in  column  five 
of  Table  6. 

6.2.2  Aluminum  1100  Alloy 

Five  1100-0  A1  specimens  were  machined  and  tested 
in  tension  at  the  University.  The  stress-strain  curves  from 
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Figure  22.  Measured  Stress-strain  Curve  for  a  1020  Steel 
Specimen  (Shot  24) . 
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Figure  23.  Averaged  Stress-strain  Behavior  for  1020  Steel  at 
Various  Strain  Rates. 
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TABLE  6 

SUMMARY  OF  TESTS  AND  RESULTS 


Shot  # 

Test 

Type 

Material 

and 

Supplier 

Strain 

Rate 

(1/s) 

Yield 
Stress 
(kbar ) 

True 

Yield 

Stress 

(kbar) 

HB016 

Tension 

1020  Steel 

370 

5. 5±0.45 

HB024 

(UD) 

730 

6.6±0.4 

!  HB026 

*« 

II 

945 

7.2±0.5 

i  HB027 

II 

If 

945 

9.1+0.47 

1  HB029 

1 

fl 

" 

975 

7. 3+0. 4 

HB032 

II 

1100  A1 

1215 

1.8±0.45 

HB033 

If 

(UD) 

1200 

2.65+0.45 

HB034 

II 

II 

890 

1.75±0.42 

HB036 

II 

M 

870 

1.96+0.5 

HB037 

It 

II 

550 

1.85+0.4 

HB040 

II 

A1  (3003) 

330 

1.95+0.2 

HB041 

If 

It 

340 

2.37+0.2 

HB043 

II 

II 

185 

1.94+0.1 

HB044 

" 

If 

555 

2.06+0.25 

HB059 

II 

II 

525 

3.12+0.5 

HB060 

11 

II 

650 

2.95+0.2 

HB064 

II 

II 

1450 

2*68j:0*  2 

HB066 

11 

II 

360 

2.010.1 

HB047 

II 

2024-T4  A1 

280 

4. 3*0. 3 

HB048 

It 

(Eglin) 

475 

4 . 5l0. 3 

HB04  9 

It 

II 

480 

4.4+0. 3 

HB050* 

11 

II 

800 

5. 1+0. 3 

HB051 

11 

OFHC  Copper 

545 

2.01+0.22 

HB052 

II 

(Eglin) 

525 

2.5+0. 3 

HB053 

It 

fl 

355 

2.45±0.2 

HB054 

It 

II 

860 

2.45+0.4 

HB045 

II 

OFHC  Copper 
(UD) 

750 

4.65+0.25 

HB056 

II 

Nickel  200 
(UD) 

770 

3.42+0.22 

HB069 

Compression 

A1  (3003) 

2900 

3. 34±0.07 

3.02 

HB070 

If 

II 

2875 

3.6+0.07 

3.2 

HB071 

1 

If 

II 

2780 

3.3+0.07 

2.98 

*  The  observed  strain  at  which  the  specimen  failed  was  0.21. 
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these  tests  indicate  that  the  yield  stress  is  slightly  sensitive 
to  the  strain  rate,  for  the  range  tested  between  550  and  1210  s 
(see  Figure  24) .  The  relatively  large  uncertainty  in  the  yield 
stress  values  evident  in  i igure  24  is  a  result  of  the  low  gain 
setting  of  the  oscilloscope  during  the  testing.  These  results 
agree  with  the  work  reported  in  references  19  and  20,  who  found 
that  high  purity  aluminum  is  highly  strain-rate  sensitive,  while 
aluminum  alloys  are  insensitive  over  the  same  range  of  strain 
rates . 

6.2.3  Aluminum  3003  Alloy 

Several  3003  aluminum  specimens  were  supplied  by 
BRL.  These  were  tested  at  the  University  in  tension  and 
compression.  Tests  were  conducted  at  various  strain-rate 
levels  using  the  split  Hopkinson  Bar  technique.  Tests  were 
also  conducted  under  quasi-static  loading  conditions  so  that 


Figure  24.  Yield  Stress  Measured  for  1100  Aluminum. 
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the  low  strain-rate  behavior  of  the  material  could  be  compared 
to  the  dynamic  behavior. 

Data  from  individual  stress-strain  tests  have 
been  averaged  in  order  to  obtain  the  characteristic  behavior 
at  various  strain  rates.  The  average  results  are  shown  in 
Figure  25.  The  strain-rate  dependency  of  the  material  is 
clearly  evident.  Considerable  increase  in  the  value  of  dynamic 
flow  stress  occurs  even  at  intermediate  strain-rate  levels 
(500  to  600  s”^) .  The  strain-rate  sensitivity  of  the  material 
seems  to  be  decreasing  with  the  increasing  strain-rate  beyond 
this  level.  Within  the  scatter  in  the  data,  the  maximiim  flow 
stress  for  a  strain-rate  of  1400/s  does  not  differ  from  the 
value  at  intermediate  strain-rate  levels.  This  may  mean  that 
variation  of  flow  stress  with  strain-rate  is  not  significant 
for  e  >  'V'500  s 


STRAIN  (%) 


Figure  25.  Average  Stress-strain  Behavior  of  3003-aluminum  at 
Various  Strain  Rates. 
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It  is  also  interesting  to  observe  at  a  strain- 

rate  level  of  1400  s  ,  the  flow  stress  reaches  maximvun  at  the 

initial  yield,  and  then  decreases.  This  can  be  interpreted 

as  analogous  to  the  strain  softening  which  occurs  at  point-of 

load  maximum  in  a  static  test.  At  high  strain  rates,  the 

elongation  is  quite  considerable  beyond  the  initial  plastic 

instability  point.  This  confirms  the  theoretical  result, 

121) 

reported  by  Hutchinson  and  Neale  ,  that  there  is  an  increase 
in  the  post-uniform  elongation  due  to  an  increase  in  strain- 
rate  sensitivity  (m)  of  the  material. 

Figure  26  is  a  schematic  representation  of  the 

data  in  Figure  25,  in  which  two  types  of  instabilities  are 

illustrated.  While  the  first-instability  point  (A)  corresponds 

to  the  load  maximum  point,  the  second  instability  point  (B) 

corresponds  to  localized  necking.  The  theoretical  results 

indicate  that  the  first  point  of  plastic  instability  occurs 

at  lower  values  of  the  uniform  strain  for  increasing  strain- 
(22  23) 

rate  levels  '  .  Figure  25  is  in  agreeraent  with  this  pre¬ 

diction.  In  the  limiting  case,  the  point  (A)  occurs  before 
any  strain  hardening  takes  place,  as  shown  by  the  1400  s  ^ 
cruve. 

The  deformation  process  beyond  the  load  maximum 

point  (A)  is  quasi-stable  and  the  amount  of  uniform  straining 

depends  on  the  strain-rate  dependency  of  the  material.  Several 

investigators  have  reported  that  ductility  in  rate  sensitive 

materials  decreases  as  strain  rate  increases.  These  predictions 

(24) 

were  experimentally  confirmed  by  Sagat  and  Tapling  for 

quasi-static  loading.  However,  under  dynamic  loading  conditions 
observed  here,  ductility  increased  with  strain-rate  level. 

This  means  that  the  inertia  of  the  deformation  process  post¬ 
phones  the  occurence  of  localized  necking  which  precipitates 
the  instability  (B)  in  Figure  26.  The  test  results  on  3003- 
aluminum  are  also  agreement  with  the  observation  of  additional 
ductility  in  copper  and  6061  aluminum  at  high  strain  rates^^^'^^1 


48 


Figure  26.  Low  Strain  Rate  and  High  Strain  Rate  Behavior, 
Illustrating  Instabilities. 

In  one  of  the  tests,  the  sample  was  dynamically 
strained  without  failure  to  an  uniform  strain  of  11.5  percent. 
The  maximum  strain  rate  in  this  test  was  200  s~^.  The  dynam¬ 
ically  pre-strained  sample  was  then  statically  loaded  to 
failure  as  shown  in  Figure  27.  The  static  portion  of  the 
stress-strain  curve  of  this  sample  apparently  follows  the 
same  static  stress-strain  curve  as  samples  without  any  dynamic 
pre-strain.  It  can  be  concluded  that  the  material  behavior 
is  independent  of  past  history  (dynamically  loaded  and  then 

unlaoded) .  Similar  behavior  was  observed  in  aluminum  6061-T4 
( 14 ) 

and  A1  6061-T6  .  These  ere  relatively  strain-rate  indepen¬ 

dent  materials,  whereas  the  present  results  were  obtained  for 
a  rate-dependent  material. 

In  contrast  to  Figure  27,  the  test  results  of  Fyfe 
( 14 ) 

and  Rajendran  for  a  strain-rate  dependent-material,  namely, 

annealed  copper  indicate  that  plastic  flow  behavior  upon 
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3.5 


30 


DYNAMICALLY  LOADED 


Fiaure  27.  Comparison  of  Dynamic,  Static,  and  Dynamic~Pre-“ 
strained  Behavior  of  3003-aluminum. 

static  reloading  on  dynamically  prestrained  thin  cylinders  does 
depend  on  the  past  dynamic  loading;  the  material  is  strain- 
rate  hardened  due  to  the  dynamic  loading.  The  static  stress- 
strain  curve  of  the  dynamically  restrained  sample  remains 
higher  than  the  static  stress-strain  curve  of  the  sample  with¬ 
out  any  prior  history. 

Three  3003  Al  specimens  were  tested  in  compression 
at  a  strain-rate  of  2850  s  The  records  from  one  of  the 
compression  tests  are  presented  in  Figure  28.  True  stress  and 
true  strain  data  were  calculated  for  these  tests  as  well. 

The  true  stress-true  strain  curves  show  an  upper  and  lower 
yield  points  and  a  relatively  constant  flow  stress  (see  Figure 
29) .  The  upper  yield,  lower  yield,  and  the  flow  stresses 
were  3.0,  1.85,  and  2.3  kbars,  respectively. 
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STRESS  (kbar) 


STRAIN 

(a) 


Figure  28.  Compression  Test  Data  for  3003  Aluminum.  (a)  Stress- 
strain  Curve,  (b)  Strain  Rate,  Strain,  and  Strain 
as  Functions  of  Time,  (c)  True  Stress-True  Strain 
Curve . 
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STRAIN  RATE  (S  ) 


TIME  (/is) 

(b) 

Figure  28.  Compression  Test  Data  for  3003  Aluminum.  (a)  Stress- 
strain  Curve,  (b)  Strain  Rate,  Strain,  and  Stress  as 
Functions  of  Time,  (c)  True  Stress-True  Strain 
Curve. 


TRUE  STRAIN 


(c) 

Figure  28.  Compression  Test  Data  for  3003  Aluminum.  (a)  Stress 
strain  Curve,  (b)  Strain  Rate,  Strain,  and  Stress  as 
Functions  of  Time,  (c)  True  Stress-True  Strain  Curve 


TRUE  STRESS  (kbar) 


TRUE  STRAIN 

(a) 

Figure  29.  True  Stress-True  Strain  Curve  for  3003  Aluminum  in 

Compression.  (a)  shot  69,  (b)  shot  70,  (c)  shot  71. 
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TRUE  STRESS  (kbar) 


(b) 

Figure  29.  True  Stress-True  Strain  Curve  for  3003  Aluminum  in 

Compression.  (a)  shot  69,  (b)  shot  70,  (c)  shot  71. 
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TRUE  STRAIN 
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Figure  29.  True  Stress-True  Strain  Curve  for  3003  Aluminum  in 
Compression.  (a)  shot  69,  (b)  shot  70,  (c)  shot  71 


6.2.4  Aluminum  2024-T4  Allo' 


Five  2024-T4  aluminum  specimens  were  supplied 
by  AFATL/DLJW  (Eglin  Air  Force  Base) .  Four  were  tested  dynam¬ 
ically  in  tension  with  the  Hopkinson  Bar  apparatus  at  strain 
rates  ranging  from  280  to  800  s  and  one  was  loaded  statically 
to  failure  at  a  0.001  s  ^  strain  rate.  Smoothed  stress-strain 
curves  of  this  material  are  presented  in  Figure  30,  and  a  true 
stress-true  strain  curve  from  the  low  rate  test  is  presented 
in  Figure  31.  These  stress-strain  curves  show  that  this 
material  is  strain  rate  sensitive.  The  yield  stress  increased 
with  strain  rate.  They  also  show  that  the  2024-T4  strain 
hardens  at  low  and  moderate  i  .rain  rates  (up  to  500  s  ■*■)  .  The 
flow  stress  becomes  constant  at  higher  strain  rates.  Shot 
HB050  was  particularly  interesting;  the  specimen  failed  during 
the  first  tensile  wave,  and  the  observed  stress  and  strain  of 
failure  were  5.1  kbar  and  0.21  respectively. 

6.2.5  OFHC  Copper 

Four  annealed  OFHC  copper  specimens  supplied  by 
AFATL/DLJW  (Eglin  Air  Force  Base)  were  tested  in  tension  at 
strain  rates  ranging  from  350  to  860  s  The  individual  data 
indicated  that  this  material  has  very  low  strain  rate  sensitivity 
for  the  range  tested.  An  average  stress-strain  curve  from  these 
tests  is  presented  in  Figure  32.  Similar  results  were  reported 
in  references  19  and  20.  A  correction  of  the  data  to  true 
stress  reveals  that  flow  stress  continuously  increases  with 
strain. 

One  OFHC  copper  specimen  was  machined  at  the 
University  and  tested  in  tension  at  a  strain  rate  of  750  s 
A  smoothed  stress-strain  curve  from  this  test  is  presented  in 
Figure  33.  The  yield  stress  of  this  specimen  was  considerably 
higher  than  the  yield  stress,  observed  at  similar  strain  rates, 
of  the  annealed  OFHC  supplied  by  AFATL/DLJW.  This  discrepancy 
is  attributed  to  the  cold  work  effects.  The  specimen  was  not 


annealed  after  machining.  This  interpretation  agrees  with 
the  effects  of  cold  work  that  was  reported  in  reference 
20. 


f 

! 


[ 

t 

I 

i 


STRESS  (kbar) 


STRAIN 


Figure  30.  Stress-strain  Behavior  of  2024-T4  Aluminum  at 
Various  Strain  Rates. 
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STRESS  (kbor) 


Stress-strain  Behavior  of  Worked  OFHC  Copper. 
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